Finite axisymmetric charged dust disks in conformastatic spacetimes 



00 

o 
o 

(N 



Guillermo A. GonzaleJB 

Escuela de Fisica, Universidad Industrial de Santander, A. A. 678, Bucaramanga, Colombia and 
Departamento de Fisica Teorica, Universidad del Pais Vasco, 48080 Bilbao, Spain 

Antonio C. Gutierrez-Pinere^j] and Paolo A. Ospin£@ 
Escuela de Fisica, Universidad Industrial de Santander, A. A. 678, Bucaramanga, Colombia 

An infinite family of axisymmetric charged dust disks of finite extension is presented. The disks are 
obtained by solving the vacuum Einstein-Maxwell equations for conformastatic spacetimes, which 
are characterized by only one metric function. In order to obtain the solutions, it is assumed that 
the metric function and the electric potential are functionally related and that the metric function 
is functionally dependent of another auxiliary function, which is taken as a solution of Laplace 
equation. The solutions for the auxiliary function are then taken as given by the infinite family of 
generalized Kalnajs disks recently obtained by Gonzalez and Reina [4J|, which is expressed in terms 
of the oblate spheroidal coordinates and represents a well behaved family of finite axisymmetric 
flat galaxy models. The so obtained relativistic thin disks have then a charge density that is equal, 
except maybe by a sign, to their mass density, in such a way that the electric and gravitational 
forces are in exact balance. The energy density of the disks is everywhere positive and well behaved, 
vanishing at the edge. Accordingly, as the disks are made of dust, their energy-momentum tensor 
it agrees with all the energy conditions. 
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I. INTRODUCTION 

The study of axially symmetric solutions of Einstein 
and Einstein-Maxwell field equations corresponding to 
disklike configurations of matter, apart from its purely 
mathematical interest, has a clear astrophysical rele- 
vance. Indeed, thin disks can be used to model accretion 
disks, galaxies in thermodynamical equilibrium and the 
superposition of a black hole and a galaxy. Disk sources 
for stationary axially symmetric spacetimes with mag- 
netic fields are also of astrophysical importance mainly 
in the study of neutron stars, white dwarfs and galaxy 
formation. Now, although normally it is considered that 
disks with electric fields do not have clear astrophysical 
importance, there exists the posibility that some galax- 
ies be positively charged [l[ , so that the study of charged 
disks may be of interest not only in the context of ex- 
act solutions. Accordingly, through the years, many at- 
tempts had been made to find exact solutions, static and 
stationary, to the Einstein and Einstein-Maxwell equa- 
tions that have as its source a relativistic thin disk. 

Exact solutions having as sources relativistic static 
thin disks were first studied by Bonnor and Sackfield 
and Morgan and Morgan [j, E|- Since then sev- 
eral classes of exact solutions correspondin g to static 
I, S 0, B S M EL IH IH and stationary [HM M, E3 
thin disks have been obtained by different authors and 
the superposition of a static or stationary thin disk with 
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a black hole has been considered in [18|, |19|, |2fJ, |21|, |22], |23j, 
HH HE Hfl- Relativistic disks embedded in an expand- 
ing FRW universe have been studied in [2tJ , perfect fluid 
disks with halos in [28| and the stability of thin disks 
models has been investigated using a first order pertur- 
bation of the energy- momentum tensor in (29[. On the 
other hand, thin disks have been discussed as sources 
for Kerr-Newman fields [3(| [Mj] , magnetostatic axisym- 
metric fields [321] a nd conformastatic and conformasta- 
tionary metrics |33l. l34l ] , while models of electrovacuum 
static counterrotating dust disks were presented in [35] ]. 
Charged perfect fluid disks were also studied in [3(| , and 
charged perfect fluid disks as sources of static and Taub- 
NUT-type spacetimes in (HHl]. 

Now, between the above mentioned works, particularly 
interesting are those that consider dust disks in comfor- 
mastatic spacetimes [33|, In this case, the charge 

density of the disks is equal to their mass density and so 
the electric and gravitational forces are in exact balance. 
This kind of equilibrium configuration has been called 
by some authors 'electrically counterpoised dust' (ECD) 
and has been studied with some detail, both in classical 
and relativistic theories [H, \M, EE S El ■ Now, as the 
matter content of the source is dust, the corresponding 
energy-momentum tensor very probably will agree with 
all the energy conditions, a fact that is particularly rele- 
vant for models of relativistic thin disks. Indeed, as can 
be see in some of the above mentioned works, there are 
many models of relativistic thin disk that do not agree 
with these conditions, and also many models that only 
fulfill these conditions partially. 

The conformastatic thin disks presented at references 
[33l [33 ] were obtained by means of the well known 'dis- 
place, cut and reflect' method in order to introduce a dis- 
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continuity at the first derivative of one otherwise smooth 
solution. The result is a solution with a singularity of the 
delta function type in all the z — hypersurface and so 
can be interpreted as an infinite thin disk. On the other 
hand, solutions that can be interpreted as thin disks of 
finite extension can be obtained if a proper coordinate 
system is introduced. A coordinate system that it adapts 
naturally to a finite source and it presents the required 
discontinuous behavior is given by the oblate spheroidal 
coordinates. Some examples of finite thin disks solutions 
expressed in these coordinates can be found at references 

alia. 

In this paper we present an infinite family of axially 
symmetric charged dust disks of finite extension. In 
order to obtain the thin disk solutions, we will solve 
the Einstein-Maxwell equations for conformastatic space- 
times assuming that the metric function and the electric 
potential are functionally related and that the metric 
function is functionally dependent of another auxiliary 
function, which is taken as a solution of Laplace equation. 
Then we take the solutions for the auxiliary potential as 
given by the infinite family of generalized Kalnajs disks 
recently obtained by Gonzalez and Reina |44|, which is 
expressed in terms of the oblate spheroidal coordinates 
and represents a well behaved family of finite axisymmet- 
ric flat galaxy models. 

Accordingly, the paper is organized as follows. First, in 
Sec. HU we present the vacuum Einstein-Maxwell equa- 
tions system and their solution for conformastatic space- 
times. We introduce the assumed functional dependences 
in order to explicitly integrate the equations system in 
such a way that the metric function and the electric po- 
tential can be expressed in terms of solutions of Laplace 
equation. The oblate spheroidal coordinates are intro- 
duced and the general solution of the Laplace equation 
expressed in these coordinates is presented. 

Next, in Sec. IIII1 we obtain the surface energy- 
momentum tensor and the surface current density of the 
relativistic thin disks. So, we first present a summary 
of the procedure to obtain the energy-momentum tensor 
and the current density by using the distributional ap- 
proach for the case of conformastatic spacetimes. The 
relation between the energy density and the charge den- 
sity is explicitly derived and the energy density is written 
in terms of the mass density of the Newtonian thin disks, 
in such a way that is explicitly shown that the relativis- 
tic thin disks will be in agreement with all the energy 
conditions. 

Then, in Sec. IIV| we restrict the general model by con- 
sidering a particular family of well behaved charged dust 
disks. We present the particular form of the solutions of 
Laplace equation that describe finite Newtonian thin disk 
with a well behaved surface mass density, and we use this 
solution to obtain the corresponding relativistic charged 
dust disks. The behavior of the energy density is then 
graphically analyzed and their main properties are de- 
scribed. Finally, in Sec. [Vj we present a brief discussion 
of our main results. 



II. EINSTEIN-MAXWELL EQUATIONS AND 
FINITE THIN DISK SOLUTIONS 



The vacuum Einstein-Maxwell equations, in geome- 
trized units such that c = G = /io = £o = L can be 
written as 



Gab — 87T Tab, 



pab 



0. 



(1) 



(2) 



with the electromagnetic energy-mementum tensor given 
by 



Tab — -r~ 
4ir 



F ac F b c — -g a bFcdl 



cd 



where 



(3) 



(4) 



is the electromagnetic field tensor and A a is the electro- 
magnetic four potential. 

Now then, for a conformastatic spacetime the line el- 
ement can be written in cylindrical coordinates x a = 
(t, <p, r, z) as 



ds 2 = 



3 2X dt 2 + e- 2X (r 2 d<p 2 + dr 2 + dz 2 ), (5) 



where the metric function A do not depends on t. So, if 
we take the electromagnetic potential as 



Aa = (-0,0,0,0), 



(6) 



where it is assumed that the electric potential <f> also is 
independent of t, the vacuum Einstein-Maxwell equations 
reduce to 



V 2 A 
X7 2 6 



e" 2A V0 • V(/», 



2 VA • V(j>, 



-2A, 



(7) 
(8) 
(9) 



where i, j = 1,2,3 and V is the usual differential operator 
in cylindrical coordinates. 

In order to solve the above equations system, first it 
is assumed that the electric potential <p is functionally 
dependent of the metric function A, (f> — cj)(X), so that 
equation ([9]) implies that 



[0'(A)] 2 = e 2 \ 
whose solution is given by 



fci, 



(10) 



(11) 



with k\ an arbitrary integration constant. With this so- 
lution, the equations system ([T])-© reduces to only one 
non-linear partial differential equation, 



V 2 A = VA • VA, 



(12) 
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for the metric function A. Then, we assume an additional 
functional dependence for the metric function A, which 
it is expressed as A = X(U), where U is an auxiliary 
function that is taken as a solution of Laplace equation. 
Thus then, equation (fT2j) it reduces to 



(13) 



x"(u) = [\'(u)Y 



whose solution is given by 



U + k 2 



(14) 



where k 2 and k% are arbitrary integration constants. 

Now, in order to have an appropriated behavior at in- 
finity, we will impose some boundary conditions on the 
solutions. So, in order to have an assymptotically flat 
spacetime, we must have that e x = 1 at infinite. Also we 
require that <p — at infinity, since we are only consider- 
ing sources of finite extension. Accordingly, we shall con- 
sider for the function U only those solutions of Laplace 
equation that behave as U — at infinite, in such a way 
that by taking k 2 = k$ we obtain the desired behavior 
for A. Finally, in order to fulfill the required behavior for 
<f>, we must take k\ = =Fl and so equation (|11[) can be 
cast as 



= ± 



U 



1 



(15) 



where the constants have been renamed as k 2 = k 3 = k. 

As we can see from (|14|) and (|15|) . the solutions of the 
Einstein-Maxwell equations system (JT])-© are expressed 
in terms of U, a solution of the Laplace equation. So, 
in order to have solutions that correspond to finite thin 
disks we must to use solutions of Laplace equation that 
properly describe Newtonian disklike sources of finite ra- 
dius. This kind of solutions can be obtained by introduc- 
ing the oblate spheroidal coordinates, which adapt in a 
natural way to the geometry of these sources. This co- 
ordinates are related with the cylindrical coordinates by 
the relation 



a a (l + £ a )(l-», 2 ), 



(16) 



(17) 



where < £ < oo and — 1 < ry < 1. So, the disk has 
coordinates £ = 0, < ?f < 1 and, on crossing the disk, 
■q changes sign but does not change in absolute value. 

Now then, the singular behavior of the coordinate r) 
implies that a polynomial in even powers of 77 is a con- 
tinuous function everywhere but has a discontinuous r\- 
derivative at the disk. Accordingly, the general solution 
of Laplace equation corresponding to an axially symmet- 
ric thin disklike source of radius a can be expressed as 

El 
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C 2n q 2n (£)P 2n {v), 



(18) 



where the P 2n (v) are the legendre polynomials of order 
In and q 2n (Q = * 2 " +1 <32n(«£), wrtn Q2n(i£) the legendre 
function of second kind of imaginary argument. The C 2n 
are constants that must be specified in order to have 
any particular solution and are directly related with the 
surface mass density of the Newtonian source considered. 
So, in a next section, we will present a particular choice 
of these constants corresponding to an infinite family of 
finite thin disks with a well behaved surface mass density, 
the family of Kalnajs generalized disks recently presented 
by Gonzalez and Reina (44[. 



III. ENERGY-MOMENTUM TENSOR AND 
CURRENT DENSITY 

As was pointed in the precedent section, the solutions 
of the Einstein-Maxwell equations corresponding to a fi- 
nite disklike source are even functions of the z coordi- 
nate. So then, they are everywhere continuous functions 
but with their first z-derivatives discontinuous at the 
disk surface. Accordingly, in order to obtain the energy- 
momentun tensor and the current density of the source, 
we will express the jump across the disk of the first z- 
dcrivatives of the metric tensor as 



bab — [gab,z] — ^9a- 



b.z 



=0+ ' 



(19) 



and the jump across the disk of the electromagnetic field 
tensor as 



[F za ] = [A a<z ] = 2A a , 



l,=0+ ' 



(20) 



where the reflection symmetry of the functions with re- 
spect to z = has been used. 

Then, by using the distributional approach [H, El, 
H3|, the Einstein-Maxwell equations yield an energy- 
momentum tensor as 



ab / 



z) + T« b [l-9{z)} + Q ab 5{z) 



and a current density as 

J a = I a 8{z), 



(21) 



(22) 



where 0(z) and 5(z) are, respectively, the Heaveside and 
Dirac distributions with support on z = 0. Here are 
the electromagnetic energy-momentum tensors as defined 
by §5§ for the z > and z < regions, respectively, 
whereas that 

16ttQ£ = b az 8 7 h - - b zz 5t + g az b z b - g zz b a h 

+ K{g zz & a b - g az S z b ) (23) 

gives the part of the energy-momentum tensor corre- 
sponding to the disk source, and 



47r/ a = [F az ] 



(24) 



is the contribution of the disk source to the current den- 
sity. Now, the "true" surface energy-momentum tensor 
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of the disk, S a b, and the "true" surface current density, 
j a , can be obtained through the relationships 



Sab 



r 



Qab S(z) ds n = e A Q ab , (25) 



I a S(z) ds n 



e- x I a , 



(26) 



where ds n = yjg zz dz is the physical measurement of 
length in the direction normal to the disk. 

For the metric ([5]), the only non-zero component of Q% 

is 



Q° = 



; 2A A„ 
2tt 



whereas that the only non-zero component of I a is 



jo = H>,* 
2tt' 



(27) 



(28) 



Thus, the only non-zero component of surface energy- 
momentum tensor Si is given by 



cO _ 
d — 



3% 

2tt 



(29) 



and the only non-zero component of the surface current 
density j a is 



e~V,, 
2tt 



(30) 



where all the quantities are evaluated at z = + . 

The surface energy-momentum tensor of the disk and 
the surface current density of the disk then can be written 
as 



where 



S ab = eV a V b 



V a = e" A (l, 0,0,0) 



(31) 



(32) 



(33) 



is the velocity vector of the matter distribution. So then, 
the energy density and the charge density of the distri- 
bution of matter are given by 



2tt 

<j>,z 

2ir 



(34) 



(35) 



respectively. Now, by using equation (jlip , the expression 
J5\i can be written as 



° = T e 



(36) 



so that the charge density of the disks is equal, except 
maybe by a sign, to their mass density. Accordingly, the 
electric and gravitational forces are in exact balance, as 
in the configurations of ECD that were mentioned at the 
introduction. 

Now then, in the context of classical general relativity, 
it is assumed that the energy-momentum tensor must 
fulfill certain requirements, which are embodied in the 
weak, strong and dominat energy conditions (5l| . Indeed, 
for the case of a dust source, all these conditions it reduce 
to the condition that the energy density be greater or 
equal to zero, e > 0. On the other hand, from equation 
([21, we have that the energy density can be written as 



{U + kf 



where 



2tt 



(37) 



(38) 



is the Newtonian mass density of a disklike source which 
gravitational potential is given by U. Accordingly, if the 
Newtonian mass density £ is everywhere no negative, 
the corresponding relativistic energy density e will be no 
negative everywhere only if we take k < 0. So then, in 
order that the energy-momentum tensor of the disks will 
agree with all the energy conditions, the Cm constants in 
(TTS")) must be properly chosen in such a way that X > 0. 
Furthermore, if we have a Newtonian potential U that 
it is negative everywhere, as is expected for a compact 
Newtonian source, then the energy density e will be non- 
singular everywhere at the disk. 



IV. A PARTICULAR FAMILY OF DISKS 

Now we shall restrict the previous general model by 
considering a particular family of disks with a well- 
behaved surface energy density. The members of the 
family are expressed in terms of particular solutions U m 
of the Laplace equation obtained by choosing properly 
the constants of the general solution l|18p . The obtained 
particular solutions U m represent the Newtonian gravi- 
tational potential of finite thin disks with mass density 
given by 



S m (r) = 



(2m + l)M 
2^2 



r 



(39) 



where M and a are, respectively, the total mass and the 
radius of the disk and we must take m > 1. For each 
value of m, the constants Cm are defined through the 
relation Q 



Co 



(2n+l)«ta+i(0)' 



(40) 



5 



TABLE I: The C 2 „ constants for m = 0, 6. 



m 


Co 


c 2 




c 4 


c 6 


C s 


Cio 


Cl2 


1 


M 

a 


M_ 

a 














2 


a 


10A/ 
7a 




3Af 
7a 










3 


a 


5 M 
3 a 




9M 
11a 


Mi 

3:-! a 








A 

-± 


M 
a 


20M 
11a 


162M 
143a 


4M 
11a 


7AI 
143a 






5 


M_ 

a 


25M 
13a 




ISA/ 
13a 


10 M 
17a 


35Af 
247a 


63Af 
4199a 




6 


M_ 

a 


2M 
a 




27A/ 
17a 


260A1" 
323a 


5M 
19a 


378 M 
7429a 


33 M 
7429a 


where 


^2„ = 


M 


r tt 1 / 2 


(4n + 1) 


(2m + 


1)! " 


(41) 






2a 


.2 2m (m-n)!r( 


m + n 


+ §)- 



for n < m, and C2n = for n > m. In TableUwe list the 
values of the nonzero Ci n for the first six members of the 
family. So, by using this constants in the general solution 
(fTSf is easy to see that the gravitational potential U m will 



be negative everywhere, as was imposed at the previous 
section. 

With the above values for the Cz n , we can easily com- 
pute the corresponding energy density of the disks by us- 
ing equation (|37|) . Now, in order to graphically illustrate 
the behavior of the different particular models, first we 
introduce dimensionless quantities through the relations 



U m (r) 
S m (f) 



M ' 

7rq 2 S] m (f) 
M 



e m {f) = nae m {f), 



(42) 

(43) 
(44) 



where f = r/a, < r < 1, and U m (f) is evaluated at 
z = + . Accordingly, the dimesionless energy density 
e(f) can be written as 



>(?) 



fcS m (f) 
[U m (f) + ~kf 



(45) 



with k = (ka)/M. 

Then, by using the above expressions and the values 
of the Ci n constants given at Table U we obtain the fol- 
lowing expressions 



J 



ei = 



£2 = 



£3 = 



£4 



3fcVl -r 2 



2[fc + f (f 2 -2)] 2 



5fc(l - f 2 ) 3 / 2 



2[fc-i|f(3f 4 -8f 2 +8)] 2 ' 



7fc(l - r 2 f' 2 



2[k + |ff (5f 6 - 18f 4 + 24f 2 - 16)] 2 ' 

9fc(l_f2)7/2 



2ffc 



3157T 



£5 



£6 = 



, , 768 -(35f 8 - 160f 6 + 288f 4 - 256f 2 + 128)] 2 

llfc(l-f 2 ) 9 / 2 

2[k + iUfa (63f 10 - 350f 8 + 800f 6 - 960f 4 + 640f 2 - 256)] 2 : 

13Jfc(l -f 2 ) 11 / 2 



2[k 



30037T 
2097152 



(231f 12 - 15127= 10 + 420(F 8 - 6400f 6 + 5760f 4 - 3072f 2 + 1024)] 2 



(46) 
(47) 
(48) 
(49) 
(50) 
(51) 



r 



for the energy densities of the first six disk models. 

We plot, in Figure [TJ the dimensionless surface energy 
density e m as a function of r for the first six disk models 
with to = 1, 2, 3, 4, 5 and 6. In each case, we plot e m (f) 



for < f < 1 with different values of the parameter k. 
We first take k = —0.25, the bottom curve in each plot, 
and then k = —0.5, —1, —2, —4 and —8. As we can see, 
in all the cases the energy density is everywhere positive 
and vanishes at the edge of the disk. However, there 
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FIG. 1: Dimensionless surface energy density e m as a function of f for the first six disk models with m = 1, 2, 3, 4, 5, 6. In each 
case, we plot e m (f) for < f < 1 with different values of the parameter k. We first take k = —0.25, the bottom curve in each 
plot, and then k — —0.5, —1, —2, —4 and —8. 



are different behaviors depending of the values of m and 
k. So, for the first two models, for m — 1 and m = 2, 
we find that for small values of |fc| the energy density 
presents a maximum near the edge of the disk, whereas 
that for higher values of \k\ the maximum occurs at the 
center of the disk. On the other hand, for m > 3, we find 
that for all the values of \k\ the maximum of the energy 



density occurs at the center of the disk. 

We can also see that as the value of m increases, the 
energy density is more concentrated at the center of the 
disks. Also, in the central part of the disks, as m increases 
the value of the energy density also increases in such a 
way that, for a given value of k, the energy density is 
grater for a greater value of m. On the other hand, at 
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r 

FIG. 2: Dimensionless energy density e m as a function of r 
with k = — 1 for the six disk models with m = 1, ... ,6. 



e m (0) 




-2 -4 -6 -8 -10 



k 

FIG. 3: Dimensionless energy density at the center of the 
disks, 2(0), as a fuction of k for the same six disk models 
previously considered. 

the border of the disks the behavior is opposite, the value 
of the energy density decreases as to increases. In order 
to see this, we plot in Figured] the dimensionless energy 
density e m as a function of f with k = —1 for the disk 
models with m = 1, ... ,6. Finally, for a given value of 
to, the value of the energy density initially increases with 
the value of |A;|, but then it reaches a maximum and then 
it decreases as \k\ increases. We depict this behavior at 
Figure [3] where we plot the dimensionless energy density 
at the center of the disks, e(0), as a fuction of k for the 
same six disk models previously considered. Now then, 
as the charge density of the disks is equal, up tu a sign, 
to their energy density, all the previous analysis apply as 
well to the behavior of the surface charge density of the 
models. 



V. CONCLUDING REMARKS 

We presented an infinite family of axisymmetric 
charged dust disks of finite extension with well behaved 
surface energy and charge densities. The disk were ob- 
tained by solving the vacuum Einstein-Maxwell equations 
system for a conformastatic spacetime. In order to ob- 
tain the solutions, a functional dependence was assumed 
between the electric potential and the metric function 
and beteween this one and an auxiliary function. The 
solutions were then expressed in terms of a solution of 
Laplace equation corresponding to a family of Newto- 
nian thin disks of finite radius, the generalized Kalnajs 
disks [4J], which describes a well behaved family of fiat 
galaxy models. 

The relativistic thin disks here presented have a charge 
density that is equal, up to a sign, to their energy density, 
and so they are examples of the commonly named 'electri- 
cally counterpoised dust' equilibrium configuration. The 
energy density of the disks is everywhere positive and 
well behaved, vanishing at the edge. Also, as the value 
of to increases, the energy density is more concentrated 
at the center of the disks, having a maximum at r = for 
all the values of \k\. However, for the first two models, 
for to = 1 and to = 2, for small values of the energy 
density presents a maximum near the edge of the disk, 
whereas that for higher values of |fe| the maximum occurs 
at the center of the disk. 

Furthermore, as the energy density of the disks is ev- 
erywhere positive and the disks are made of dust, all the 
models are in a complete agreement with all the energy 
conditions, a fact of particular relevance in the study of 
relativistic thin disks models. Indeed, as was mentioned 
at the introduction, many of the relativistic thin disks 
models that had been studied in the literature do not 
fully agrees with these conditions. 

Now then, as we can see from the equations at Sec. [Til 
the procedure here presented can be applied not only for 
axisymmetric conformastatic spacetimes but can also be 
used to obtain non-axisymmetric solutions of the vacuum 
Einstein-Maxwell equations. Accordingly, the thin disks 
models here presented can be generalized by considering 
for the auxiliary function U solutions of Laplace equa- 
tions without the impossed axial symmetry. So, we are 
now working in this direction and the results will be pre- 
sented in a next paper. Analogously, the generalization 
to conformastationary spacetimes with magnetic fields is 
in consideration. 
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